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ABSTRACT 



Quantum interference phenomena manifests itself in several ways in the polarized solar spec- 
trum formed due to coherent scattering processes. One such effect arises due to interference 
between the fine structure (J) states giving rise to multiplets. Another effect is that which arises 
due to interference between the hyperfine structure (F) states. We extend the redistribution ma- 
trix derived for the J-state interference to the case of F-state interference. We then incorporate 
it into the polarized radiative transfer equation and solve it for isothermal constant property 
slab atmospheres. The relevant transfer equation is solved using a polarized approximate lambda 
iteration (PALI) technique based on operator perturbation. An alternative method derived from 
the Neumann series expansion is also proposed and is found to be relatively more efficient than 
the PALI method. The effects of PRD and the F-state interference on the shapes of the linearly 
polarized Stokes profiles are discussed. The emergent Stokes profiles are computed for hypothet- 
ical line transitions arising due to hyperfine structure splitting (HFS) of the upper J = 3/2 and 
lower J = 1/2 levels of a two- level atom model with nuclear spin I s = 3/2. We confine our 
attention to the non-magnetic scattering in the collisionless regime. 



Subject headings: line: formation 
Sun: atmosphere 

1. Introduction 



methods: numerical - polarization - radiative transfer - scattering 



The linearly polarized solar spectrum is pro- 
duced by coherent scattering processes taking 
place in the solar atmosphere. This so called 
'second solar spectrum' is highly structured and 
reveals various physical processes responsible to 
generate the polarized signals in the spectrum. 
Quantum interference is one such physical process 
whose importance has been highl ighted in the 
second sola r spectrum studies (see lstenfld[i~980l 



1994 119971) . The coherent superposition of the 
fine structure states leads to the J-state interfer- 
ence, whereas the F-state interference arises due 
to superposition of the hyperfine structure states 



(see Figured]). The J-state interference theory for 
the case o f frequency coherent scatt ering was de- 
veloped by [Stenfk] ([l983[l99l[l993). This theory 
was extended to include partial freq uency redistri- 
bution (PRD) in line scattering, by ISmitha et al.l 
(|2011al hereafter called PI). The J-state PRD 
matrix derived in P I is used in the polar ized line 
transfer equation in lSmitha et al. I (l2011bl . hcrafter 
P2). An alternative scattering theory of J-state 
interference based on metalevel approa c h wa s 
developed by Landi Degrinnocenti et ah! (1997), 
which also includes the F-state interference ef- 
fects. All the papers mentioned so far are appli- 
cable to the case of colisionless regime. 
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Second solar spectrum contains several lines 
which have signatures of F-state interference. Ex- 
amples of these lines are Na I D2 at 5890 A, Ba 11 
D 2 at 4554 A, Mn 1 8741 A, Sc 11 4247 A etc. In 
this paper we are concerned with the line forma- 
tion studies involving F-state interference process 
and PRD. The F-state redistribution matrix de- 
rived in this paper can be used for modeling the 
non-magnetic quiet region observations of hyper- 
fine structure splitting (HFS) in the lines men- 
tioned above. 

The F-state interference theory applicable 
to the fre quency c oheren t scattering was de- 
veloped bv IStenflol (jl997l ). This th e ory, a long 



with PRD, was a 



pplied bv lFluri et all (|2003l ) and 



Holzreuter et al. ( 20051) in the polarized line trans- 



fcr computations. In lLandi Degl'Innocenti fc Landolfi 



(|2004l ) the theory of F-state interference was de- 
veloped under the approximation of complete fre- 
quency redistribution (CRD). The theory of F- 
state interference in a magnetic field for multi- 
term atoms in the collisionless regime and un- 
der the approximation of C RD is presented in 
Casini fe Manso Sainz] (|2005h . 

In the present paper we extend the J-state in- 
terference theory presented in PI to the case of F- 
state interference. The F-state redistribution ma- 
trix is derived here for the non-magnetic case and 
in the collisionless regime. The reason for consid- 
ering the non-magnetic case in this paper, is that 
the formulation of PI was confined to the linear 
Zeeman regime of field strengths (the spacing be- 
tween the Zeeman m-states being smaller than the 
spacing between the fine structure states). In the 
present context, the hyperfine splitting becomes 
comparable to the Zeeman splitting even for weak 
magnetic fields, and we quickly enter the Paschen- 
Back regime of field strengths (level crossing of the 
m-states belonging to different F-states) , in which 
the formulation presented in PI is not valid. Since 
the Paschen-Back effect is outside the scope of the 
present paper, our treatment here is limited to 
the non-magnetic case, but the extension to the 
Paschen-Back regime is planned to be pursued in 
future work. We further assume that the lower 
level is unpolarized and infinitely sharp. While 
this assumption is made for the sake of mathe- 
matical simplicity, it is physically justified for the 
long-lived ground states, which are correspond- 
ingly more vulnerable to collisional depolarization. 



Following P2, this PRD matrix is incorporated 
into the polarized line transfer equation, and 
solved using an operator perturbation method. 
We also propose a new method to solve the F- 
state interference problem. It is called the scatter- 
ing expansion met h od (SEM) and is d escrib ed in 
iFrisch et al.l (|2009l ) : ISampoorna et all (|201lh . Re- 
cent ly it has been applied to a variety of prob lems 
(see lSowmva et al.ll2012tlSupriva et al.ll2012l ). We 
compare the operator perturbation method and 
the SEM by applying them to the problem at 
hand. 

In Section [2] we derive the PRD matrix for F- 
state interference and incorporate it into the line 
transfer equation. In Section[3]we describe the nu- 
merical methods used to solve the transfer equa- 
tion. Results are presented in Section |U Section [5] 
is devoted to the concluding remarks. 

2. Basic equations 



2.1. The redistribution matrix 

In this section we present the redistribution ma- 
trix for the F-state interference, derived starting 
from the Kramers-Heisenberg formula. We restrict 
our attention to the case of non-magnetic collision- 
less regime. 

The redistribution matrix for the F-state inter- 
ference can be derived through a straight forward 
replacement of quantum numbers, in the J-state 
interference redistribution matrix d erived in PI . 
The replacements are as follows (see Stenflol 1997 ; 



Landi Degl'Innocenti fc Landolfill2004l ): 



L -)• J; J -> F; S -> I s , 



(1) 



where L, J and S represent the orbital, electronic, 
and spin quantum numbers of a given state. F is 
the total angular momentum and I s the nuclear 



J h ,L 



J a, I, 



Fig. 1. — Level diagram representing the HFS in a two- 
level atom model. 
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spin of the atom under consideration. The ex- 
pression for the F-state interference redistribution 
matrix expressed in terms of irreducible spherical 
tensors can be written as 



sence of a magnetic held is given by 



3(2J fc + 1) 



KF a F f F b F t 



(-1) 



21, 

Ff-Fa 



1 

cos /3 Fh ,F b e 



X )F a F s + Kf¥ b ,F„, )F a F f 

x(2F a + l)(2F / 



Jo 

F b 

Ja 
F,y 

1 



■h 
Ff 

Jb 
Ff 
1 

F b 



1 

I. 
1 

Is 

K 
F, 



1)(2F 6 - 

Ja 

F b 

Ja 

F v 
1 

F b > 



l)(2F b , 

J b 1 
F n L 



Jb 
F a 

1 

F b 



1 

I, 

K 

Ft 



xT K (i,n)T K U,n') 



(2) 



In the above expression the angle PF b ,F b is defined 
as 

tan/3 n , Ft = (3) 

7 

where fad p b , F b represent the energy differences be- 
tween the F b i and F b states in the absence of a 
magnetic field. 7 is the damping parameter of 
the upper state. The lower levels are assumed to 
be infinitely sharp and unpolarized. The h and 
/ functions are the auxiliary quantities defined 
in the same way as Equations (14) - (15) of PI, 
but with the replacements given in Equation (fTJ. 
Tq are the irreducible tensors for polar imetry 
introduced by Landi Degl'Innocenti ( 19841) . For 
the non-magnetic case presented in this paper 
Q = 0. The indices % and j refer to the Stokes 
parameters = 0, 1, 2, 3) with K = 0, 1, 2 and 
— K < Q < +K. The directions of the incoming 
and scattered rays are given by n' and n respec- 
tively, n — (9, ip) with 9 being the colatitude and 
tp being the azimuth of the outgoing ray. x' and 
x are the incoming and scattered frequencies in 
Doppler width units. 

2.2. The polarized line transfer equation 

The one dimensional radiative transfer equation 
for solving the line formation problems with PRD 
and F-state interference in scattering in the ab- 



dI(T,x,fj.) 
1 dr 



t>HFs(%)+r)[I(T, x, h)-S(t, X, (J,)], 

(4) 

where I — (I, Q) T is the Stokes vector and S = 
(Si,Sq) t is the Stokes source vector. Equa- 
tion Q is valid for the case of a two-level atom 
with an infinitely sharp and unpolarized ground 
level. /1 = cos 9 represents the line of sight, r is 
the ratio of continuum to the frequency-integrated 
line absorption coefficient. The positive Stokes Q 
represents electric vector vibrations perpendicular 
to the solar limb, r is the line optical depth de- 
fined by dr = —k^dz, where hi, is the frequency- 
integrated line absorption coefficient defined for 
a two- level atom with HFS. If r]i is the line ab- 
sorption coefficient then for the standard two-level 
atom without HFS, r\i = hi$(x) where <j)(x) is the 
Voigt profile function. In the presence of HFS, i}i 
is given by (see Equation (7) of P2) 



f\h{v) 



hi 



(21. + 1) 

Jb Ja 1 



- ]T(2F Q + l)(2F fc + l) 



F a F b 
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' < F„ F b I s ' ^ yF >>F a - v), (5) 



where 



hi 



hp 



J b Ja 



47T 



N(J a )B(J a -»• J b ), 



(6) 



is the frequency-integrated absorption coefficient 
for all the F-states. Thus ^hfs(^) is the weighted 
sum of the individual Voigt profiles 4>{vF b F a — v) 
representing each of the F a — > F b absorption. 

For the particular case of J a = 1/2 — > J b = 
3/2 — > Jf = 1/2 transition with I s = 3/2, ^>HFs(a ; ) 
takes the form 



4>hfs(x) 



2 5 
— (j){v i-u)+ 32^(^1 1 - v) 



32 



1 

32' 
14 



V12 — V) 



(7) 



We have verified that if the F-states are very 
closely spaced, then a single profile function 
4>(yj b j a — v), corresponding to the J a — > J b 
transition, can be used instead of ^H Fs(ff) (see 
Landi Degl'Innocenti fc Lan dolfi 2004). 
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The total source vector S in Equation ((4]) is 
given by 

S(t, x, /i) = -r^— , (8 

where the unpolarized continuum source vector 
S c = BU, with B being the Planck function and 
U = (1, 0) T . The line source vector for a two-level 
atom with HFS is given by 



Si(T,x,fi) = eBU 

r+l 



1 



-Hoc 



daf 



</Ws(a;) 

|^^R(x, M ; a/, *',//)• (9) 



Here e = rj/(rj + Tr) is the photon destruction 
probability per scattering also known as the ther- 
malization parameter, with Tj and Tji being the 
inelastic and radiative de-excitation rates of the 
upper state Ff,. To a first approximation these 
rates are assumed to be the same for all the in- 
states. R(x, fi; x 1 , fjf) is the redistribution matrix 
defined in Equation ([2]) but integrated over the az- 
imuths ip' of the incoming radiation. Such a sim- 
plification is possible due to the azimuthal sym- 
metry of the problem. This redistribution matrix 
can be rewritten as 



I^ij (x^ ^i, X , f-l 



K 



x7^(i,M)7^(jV)-(10) 

The redistribution function components 1Z (x, x') 
are given by 
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For simplicity, we use the angle-averaged ver- 
sions of the auxiliary functions (hp b Fb ,)F a F f and 

(f¥ b ,F b ,)F a F r 



2.3. Decomposition of the Stokes vectors 
into the reduced basis 

Decomposition of the Stokes source vector S in 
the reduced basis makes it independent of 9. The 
decomposition of S defined in Equation flSJ can be 
carried out in a way similar to the one presented 
in Section 2.1 of P2. Hence we do not repeat them 
here. The transfer equation for the reduced Stokes 
vector X can be written as 

9X(t : x : [i) 

M 7T- 2 — = (<£hfs(z) + r)[T(T, x, fx) - S(t, x)]. 

OT 

(12) 

The corresponding irreducible total and line 
source vectors are given by 

C / s (j)nFs{x)Si{T 1 x)+rg(r) 
S{t,x) = -— , (13) 



4>ufs{x) + r 



and 



/+oo 
-OO 



n(x,x') 

(t>HFs(x) 



J(t, x')dx' . 



(14) 

Here 1Z(x, x') is a (2 x 2) diagonal matrix with el- 
ements 1Z =diag (TZ Q ,1Z 2 ), where 1Z K are defined 
in Equation (|TT|) . Q(t) — (B,0) T is the primary 
source vector, and JT(t, a;) is the mean intensity 
defined in Equation (22) of P2. 

3. The numerical methods 

Here we describe two numerical techniques to 
solve the reduced form of the transfer equation. 
We compare their performance on some bench- 
mark problems. 

3.1. Operator perturbation method 

The solution of the polarized line transfer equa- 
tion defined in Equation (|12[) using the polarized 
approximate lambda iteration (PALI) method is 
described in Sections 3.1 and 3.2 of P2. The same 
equations hold good for the present problem also. 
Hence we do not repeat those equations here. The 
only difference is that the redistribution matrix 
for J-state interference is now to be replaced by 
the corresponding redistribution matrix for the in- 
state interference presented in this paper. Also the 
profile function is to be replaced with (j^uFsix). 
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3.2. Scattering expansion method 



ponent is given by 



In recent years a new method based on Neu- 
mann series expansion of th e polarized source vec- 
tor has been developed (see lFrisch et al. I I20091) . It 
is applied to a variety of astrophysical problems. 
Here we describe the application of this method to 
the problem at hand. 

In this method, the reduced line source vector 
defined in Equation (|14[) is rewritten in the com- 
ponent form for the non-magnetic case as 



S£(t,x) = G(t)Woo + 



- 1 df/_ 

-i 2 



+ 00 



' dx ' \ ( T? E***V)ffV,sV). (is) 

are the components of the Rayleigh phase 
matrix in t he reduced basis (see Appendix A of 
Frischll2007l ). We first consider the component Sq. 
Expanding the summation over K' on the right- 
hand side of Equation f(T5|) we obtain 



d 4- r dx 
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The degree of linear polarization arising due to 
Rayleigh scattering is small because of small de- 
gree of anisotropy prevailing in the solar atmo- 
sphere. Hence the effect of linear polarization 
on Stokes I can be neglected to a good approx- 
imation. Neglecting the contribution from Jg, in 
Equation (fTS|) we get 



S°(t,x)^G(t) + 



+ 1 A„l 



dfi' 



+oo 



dx 



,n°(x,x') 

</>hfs(z) 



*8V)Jofo*V),(17) 



where Sq denotes the approximate value of Sq. 
It is the solution of a non-LTE unpolarized ra- 
diative transfer equation and is computed using 
the Frequency-by-Freq uency (FBF) technique of 
Paletou fc Auen (Il995h . 

The polarization is computed from the higher 
order terms in the series expansion. The S 2 com- 



1 dfi' f +0 ° , ,K 2 {x,x') 
—— I dx 



^HFS(^) 



x^ 2 ( M ')/o 2 (r^',/i') 



(18) 



Retaining only the contribution from I® on the 
right-hand side of Equation (IT8l) we obtain the 
single scattering approximation to the polarized 
component of the source vector as 



j TZ 2 (x,x') 

X^V)^, *',//)• (19) 

The superscript (1) denotes single (first) scatter- 
ing. This solution serves as a starting point for 
the computations of higher order scattering terms. 
Thus the iterative sequence of SEM can be repre- 
sented by 

[S 2 (t,x)]M^[S$(t,x)]W 

- 1 dp,' r + °° , ,v?(x,x') 

—— j dx — t— -- 

L 2 J_ 00 </>HFSW 

x9™(n')[T>(T,x',ij!)]< n - 1 \ (20) 

where the superscript (n) denotes the n th scat- 
tering. The iterative cycle is continued until the 
required convergence criteria is met. 

In the following we compare the performance 
of these two numerical methods by plotting the 
maximum relative correction defined as 



,(*») _ 



where 



5Sf'(T,x,fj,)\ 



-s 



(n) _ 

C l ~ maX t,x,ix} -(„) 



and 



» 



max 



P( n \x,ii) - P (n - X \x,ii) 
P^-^ix^) 



(21) 



(22) 



(23) 

as a function of the iteration number as shown in 
Figure [2] In the above equations P = [Q/I] is the 



5 



degree of linear polarization and S/ = | [Sj 

Figure [5] is computed for a test problem de- 
fined by the model parameters (T, a, e, r, F) = 
(2 x 10 10 , 2 x 1CT 3 , 10~ 4 , 0, 1 ) where T is the opti- 
cal thickness of the self emitting slab and a is the 
damping parameter of the upper level J&. From 
the figure one can clearly see that the conver- 
gence rate of the SEM is larger by several fac- 
tors compared to the PALI method. The rea- 
son for the PALI method being slow is that the 
source function corrections are computed itera- 
tively from an approximate initial guess and then 
the approximate lambda operator is perturbed un- 
til the source function corrections fall below a con- 
vergence criteria. On the other hand, the ini- 
tial guess in the SEM for polarized line forma- 
tion is the single scattered solution itself (which 
already contains the physical characteristics of the 
scattering mechanism under consideration). For 
this reason SEM takes just a few iterations to 
converge to the same level of accuracy as the 
PALI method. Further, SEM is easy to imple- 
ment for problems of any physical and/or numer- 
ical complexity. This makes the SEM a method 
of choice. For a d etailed comparison of P ALI 
and SEM we refer to Sampoorna et al. ( 201 if ) and 
Supriva et ail (|2012t) . The simple Lambda itera- 



tion for polarization and the SEM are essentially 
similar. In the lambda iteration, a source vec- 
tor correction is computed at each iteration, and 
the current source vector is updated until con- 
vergence is reached. In the SEM, each iteration 
can be seen as contributing a higher order scat- 
tering term to the series expansion of polarized 
component of the source vector. This component 
is updated by adding successively higher order 
terms in the scattering expansion of the source 
vector. These points are clearly explained resp ec- 
tively in Truiillo Bueno fc Manso Sainzl ( 19991 see 
the discussion follo wing their Equation (28)), and 
iFrisch et~aL (2009, see the discussion following 
their Equation (36)). 

4. Results and discussion 

In this section we present the results computed 
for a standard two-level atom model with F-state 
interference using the PRD matrix presented in 
this paper. Isothermal constant property media 



characterized by (T, a, e, r, B) are used. The slabs 
are assumed to be self-emitting. 

The results are presented for transitions cen- 
tered at hypothetical wavelengths arising due to 
HFS of the J b = 3/2 and J a = 1/2 levels of a two- 
level atom with nuclear spin I s — 3/2. Due to the 
hyperfinc interactions the upper J-state splits into 
four F-states with F, = 0,1,2,3, and the lower 
J-state splits into F a = 1,2. Owing to the se- 
lection rule AF = 0,±1, these F-states produce 
six radiative transitions involving them (see Ta- 
ble Q]). For simplicity the Doppler width of all 
the lines is taken to be constant at AA_d = 25 
mA. In the transfer computations, a grid reso- 
lution of (N d ,N x ,N^) = (5, 417, 5) is generally 
used, where Nd is the number of depth points per 
decade in the logarithmically spaced t— grid. The 
first depth point is taken as r min = 10~ 2 . N x is 
the total number of frequency points covering the 
full line profile. is the number of co-latitudes 
0(n), taken as the 5 points of a Gauss-Legendre 
quadrature. 

4.1. F-state interference effects in the case 
of single scattering 

In this section we study the behavior of the 
F-state interference PRD matrix derived in Sec- 
tion 12.11 by computing the scattered profiles in a 
single scattering event. The results in Figure[3]are 
computed for a 90° single scattering event. This 
is done by giving as input an unpolarized beam of 
light incident on the scattering atom at (if = — 1 
and observing the scattered ray at fj, = in the 
scattering plane (see PI for details on comput- 
ing polarization profiles in a 90° single scattered 
event). The dashed line in Figure[3]is computed by 
ignoring the interference effects, whereas the solid 
line is computed by taking account of the interfer- 
ence effects between the F-states. Profile similar 



to the solid line can also be seen in Fluri et al 



(|2003l) and iHolzreuter et all (|2005l) where plots 
of the wavelength dependent polarizability factor 
W2(X) are shown. In the single scattering case, 
the profiles of the W2(A) and the Q(X)/I(X) are 
similar in shape but differ only in magnitude (see 
below) . 
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4- 1.1. Principle of spectroscopic stability for F- 
state interference 

It is well known that the principle of spectro- 
scopic stability provides a useful tool to check any 
theory of quantum interference. This was first 
discussed in the context of scatte r ing po larization 
and ap pl ied in detail in IStenflol (Il994l) (see also 



1997; Landi Degl'Innocenti fc Landolfi 



Stenflol 

2004). In this paper, we apply it to the case 
of F-state interference arising due to the nuclear 
spin I s . According to the principle of spectro- 
scopic stability, in the limit of vanishing HFS in 
a two-level atom, the theory of F— state interfer- 
ence should reduce to the standard two-level atom 
theory without HFS. This can be verified by com- 
puting the polarizability factor Wi and in turn the 
fractional polarization Q /I in the limit of vanish- 
ing F-states. The value of W 2 in this asymptotic 
limit (which can be obtained by ne glecting, t he I s ) 
can be computed as described in IStenflol (|1997h 
with 



(H^asym — 



1 1 2 

Jb Jb Ja 



1 1 2 

Jb Jb Jf 



1 1 

Jb Jb Ja 



1 1 

Jb Jb Jf 



(24) 



For the particular case of J a = 1 /2 — > = 3/2 — !> 
Jf = 1/2 scattering transition, (W2)asym = 0.5. 
Hence W 2 (A) is expected to appr oach 0.5 in th e 
very far wings (see Figure 2 of IStenflol fl997h . 
In the 90° single scattering case, the Q/I and 
the W~2(A) are related through the fo rmula (see 
Landi Degl'Innocenti fc Landolfi! [200l 



Q(X)/I(X) 



3VF 2 (A) 
4- W 2 (A)' 



(25) 



The above formula gives in the far wings a value 
of Q/I = 0.428 for (VF 2 ) asym = 0.5. 

From Figure |31 we can see that the solid curve 
reaches an asymptotic value of 42.8% as demanded 
by the principle of spectroscopic stability, whereas 
the dashed line reaches about 10% in the far wings, 
thereby violating the principle of spectroscopic 
stability. These arguments show that in the for- 
mulation of the redistribution matrix, the inclu- 
sion of interference effects between the F-states is 
essential. 



4.2. Effects of F-state interference in mul- 
tiply scattered Stokes profiles 

In this section we present the results obtained 
by solving the transfer equation including the in- 
state interference. In the particular case of opti- 
cally thin slabs, it can be shown, by choosing ap- 
propriate geometric arrangement, that the multi- 
ply scattered solution approaches single scattered 
solution thus proving that we have correctly incor- 
porated the F-state redistribution matrix in the 
line transfer code. See P2 for more details regard- 
ing single scattering in a thin atmospheric slab. 

When the optical thickness of the medium is 
large, multiple scattering effects come into play. 
Figure [4] shows one such example, where the emer- 
gent Stokes profiles are computed for different op- 
tical thicknesses. The dashed line in this figure 
is computed by neglecting HFS. This is the stan- 
dard two-level atom case which results in a single 
radiative transition. The dotted line is computed 
with HFS but without interference between the 
F-states. In this case the six radiative transitions 
arising due to HFS are treated independently. The 
solid line is computed taking account of the F- 
state interference. This comprises of six interfering 
radiative transitions between the F-states. The 
three line types in this figure are quite similar to 
each other in shape but differ prominently in am- 
plitude. 

For T = 2, the atmospheric slab is effectively 
thin and the Q/I profiles for both solid and dot- 
ted lines have a structuring within the line core 
which is different from that of the dashed line. 
This is due to the HFS of the given J-level. As 
the optical thickness increases, such a structuring 
gets smoothened out and the shape (not the ampli- 
tude) of the solid and dotted line profiles resemble 
more closely with the dashed line profiles. 

In the case of effectively thick atmospheric slabs 
(T > 2), two peaks are seen on either side of 
the line center arising due to PRD effects and are 
known as PRD wing peaks. In the line core, the 
solid and dotted lines nearly coincide whereas the 
dashed line differs from these two. This shows 
that the depolarization in the line core is purely 
due to HFS, irrespective of the interference effects 
between the F-states being included. In the wings, 
the solid line and the dashed line coincide whereas 
the dotted line differs significantly. Upon compar- 
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ing the solid and dotted lines, it is evident that 
the interference effects show up in the line wing 
PRD peaks like in the case of J-state interference. 
However the J-state interference effects are seen 
even beyond the PRD wing peaks unlike the case 
of F-state interference. When F-state interference 
is taken into account the Q /I in the wings reaches 
the value of the single line case as expected from 
the principle of spectroscopic stability (see Sec- 
tion But when interference is neglected, 
the dotted and dashed lines differ considerably in 
the wings which can be seen as a violation of the 
principle of spectroscopic stability. Thus the prin- 
ciple of spectroscopic stability serves as a powerful 
tool to check the correctness of our formulation not 
only in the case of single scattering but also in the 
radiative transfer computations. 

Though such significant signatures of HFS and 
F-state interference are seen in Q/I, the intensity 
I remains unaffected by these effects. 

4.3. Comparison with wavelength depen- 
dent polarizability theory of Stenflo 

In this section we compare our redistribution 
matrix approach and the wavelength dependent 
polarizability H^(A) theo ry for the case of F-state 
in terference p r esent ed in Stenfiol (Il997) and use d 
in lFluri et all (|2003h and Irlolzreuter et al.l (|2005h . 
The comparison is shown in Figure [5j The dotted 
lines show the profiles computed using the exact 
PRD F-state interference theory presented in Sec- 
tion[2] This is our redistribution matrix approach. 
The dashed lines show the profiles computed us- 
ing the Wi (A) approach. The values of the W2 (A) 
are calculated from Equation (|25|) using the (Q/I) 
plotted in Figure [3] (solid line). To use the W2(A) 
in radiative transfer computations we replace the 
redistribution matrix 1Z (x, x') in Equation 

by 



W^(A)[F n - A (3/2^1/2)], 



(26) 



where R u ~ A (Jb — > Jf) is the an gle-averaged fre - 
quency redistribution function of iHummerl (|l962t ) 
for a line centered at Xj b j, corresponding to the 
•h Jf transition. For the hypothetical case 
under study, we have assumed the F-states to 
be very closely spaced. Under such an assump- 
tion, a single redistribution function computed for 
the J = 3/2 — > 1/2 transition can be used to 
represent all the F-state transitions. However if 



the F-states are far apart then the redistribution 
function needs to be computed for each of the 
F) — > Ff transition. In such a case, the redistri- 
bution matrix TZ (x, x') takes the following form 
in the ^(A) approach: 



w K (x) J2 l R 



II- A 



F 



J) 



(27) 



The polarizability factor Wo (A) = 1, and H^(A) 
is the wavelength-dependent W2 factor calculated 
from Equation (|25|) . For the closely spaced F- 
states a common absorption profile function <j>(x) 
corresponding to the J a = 1/2 — > Jb = 3/2 tran- 
sition is used. But in the case of widely spaced 
F-states, the <fi(x) has to be taken as the sum of 
all the individual F a — > Fj absorption profile func- 
tions. As seen from Figure [5] both the redistribu- 
tion matrix approach and the W2 (A) approach give 
nearly the same results. 

5. Conclusions 

In the present paper we have extended the J- 
state interference formulation discussed in PI and 
P2 to the case of F-state interference. The treat- 
ment is restricted to the collisionless and non- 
magnetic reg ime. The de composition technique 
presented in Frischl ( 2007 ) is applied to the in- 
state interference problem. It helps to incorpo- 
rate the F-state interference redistribution matrix 
into the reduced form of the line radiative trans- 
fer equation. The transfer equation is solved us- 
ing the traditional PALI and the scattering expan- 
sion method by suitably adapting them to han- 
dle the F-state interference problem. The SEM 
is found to be more efficient and faster than the 
PALI method. 

The Stokes profiles computed by taking account 
of HFS are similar to the profiles of a single line 
arising from a two-level atom model without HFS. 
The HFS causes a depolarization oiQ/I in the line 
core irrespective of whether the F-state interfer- 
ence is taken into account or not. Like the J- 
state interference, the F-state interference affects 
mainly the line wing PRD peaks. We also show 
that when interference effects are neglected, the 
principle of spectroscopic stability is violated in 
both single scattered and multiple scattered pro- 
files. Using the fractional polarization Q/I in the 
90° single scattering case, we can numerically esti- 
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mate the wavelength dependent polarizability fac- 
tor W^(A). The W2(X) so computed can then be 
used in the transfer equation to compare with our 
exact redistribution matrix approach. The two 
approaches are found to give identical emergent 
Stokes profiles. 

REFERENCES 

Casini, R., & Manso Sainz, R. 2005 ApJ, 624, 1025 

Fluri, D. M., Holzreuter, R., Klemcnt, J., & Sten- 
flo, J. O. 2003, in ASP Conf. Ser. 307, So- 
lar Polarization III, ed. J. Trujillo-Bueno & J. 
Sanchez Almeida (San Francisco, CA: ASP), 
263 

Frisch, H. 2007, A&A, 476, 665 

Frisch, H., Anusha, L. S., Sampoorna, M. & Na- 
gendra, K. N. 2009, A&A, 501, 335 

Holzreuter, R., Fluri, D. M., & Stcnflo J.O. 2005, 
A&A, 434, 713 

Hummer, D. G. 1962, MNRAS, 125, 21 

Landi Degl'Innocenti E. 1984, Sol. Phys., 91, 1 

Landi Degl'Innocenti, E., Landi Degl'Innocenti, 
M., & Landolfi, M. 1997, in Proc. Forum 
THEMIS, Science with THEMIS, ed. N. Mcin & 
S. Sahal-Brechot (Paris: Obs. Paris-Meudon), 
59 

Landi Degl'Innocenti, E., & Landolfi, M. 2004, 
Polarization in Spectral Lines (Dordrecht: 
Kluwer) 

Paletou, F., & Auer, L. H. 1995, A&A, 297, 771 

Sampoorna, M., Nagendra, K. N., & Frisch, H. 

2011, A&A, 527, A89 

Smitha, H. N., Sampoorna, M., Nagendra, K. N., 
& Stenflo, J. O. 2011a, ApJ, 733, 4 (PI) 

Smitha, H. N., Nagendra, K. N., Sampoorna, M., 
& Stenflo, J. O. 2011b, A&A, 535, 35 (P2) 

Sowmya, K., Nagendra, K. N., & Sampoorna, M., 

2012, MNRAS, 423, 2949 

Stenflo, J. O. 1980, A&A, 84, 68 



Stenflo, J. O. 1994, Solar Magnetic Fields: 
Polarized Radiation Diagnostics (Dordrecht: 
Kluwer) 

Stenflo, J. O. 1997, A&A, 324, 344 

Supriya, H. D., Nagendra, K. N., Sampoorna, M., 
& Ravindra, B. 2012, MNRAS, 425, 527 

Trujillo Bueno, J. & Manso Sainz, R. 1999, ApJ, 
516, 436 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



9 



Table 1: Wavelengths (A) of F-state transitions for a hypothetical atomic system 
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Fig. 2. — Comparison of PALI (solid line) and scattering 
expansion method (dotted line). The model parameters are 
given in the text. A convergence criteria of 10 -8 is used. 
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Fig. 3. — The profiles of the intensity / and the fractional 
polarization Q/I, plotted for a hypothetical line system 
with hypcrfine structure splitting. Solid line represents the 
Q/I with F- state interference and dashed line represents 
Q/I without _F-state interference. Single 90° scattering is 
assumed at the extreme limb (fi = 0). The model parame- 
ters arc a = 0.002, the Dopplcr width AA D = 0.025 A. 
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Fig. 5. — Comparison between the redistribution matrix 
theory (dotted line) and wavelength dependent polarizabil- 
ity factor VK2(A) theory of Stenflo (dashed line). The opti- 
cal thickness of the atmospheric slab is T = 2 X 10 4 . The 
other model parameters are the same as in Figure [4] 



Fig. 4. — Comparison between the multiply scattered 
emergent Stokes profiles computed for different atomic sys- 
tems as indicated in panel 2. The model parameters are 
(a,€,r,B) = (2 X 10" 3 , 10" 4 , 0, 1). The line of sight is 
given by fi = 0.047. The wavelength positions of the six 
components are given in Table 1. The spacing between the 
hyperfine structure components is taken to be the same as 
those corresponding to the Nal D2 line. 



11 



